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Functional renormalization group study of the Nambu-Jona-Lasinio model at finite 
temperature and density in an external magnetic field 
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institute for Theoretical Physics, Kanazawa University, Kanazawa 920-1192, Japan 

In this study, we investigate the Nambu-Jona-Lasinio (NJL) model at finite temperature and finite 
density in an external magnetic field using the functional renormalization group. We investigate the 
dependence of the position of the ultraviolet fixed point (UVFP) of the four-Fermi coupling constant 
on the temperature, density, and external magnetic field, and we obtain the chiral phase structure. 

The UVFP at low temperature and finite chemical potential oscillates in a small external magnetic 
field, which can be interpreted as the de Haas-van Alphen effect. We also obtain phase diagrams 
with complex structures, where the phase boundary moves back and forth as the external magnetic 
field increases in the low temperature and high density region. 


I. INTRODUCTION 


The phase diagram and equation of state for quantum 
chromodynamics (QCD) have been studied as important 
subjects in the area of elementary particle physics. In 
particular, QCD matter in an external magnetic field 
has attracted much attention for several reasons. It 
is known that neutron stars, which are high baryonic 
density stars, generate strong external magnetic fields 
{B ~ 10 14 Gauss) 0- In heavy-ion collision experiments, 
quite strong external magnetic fields (B ~ 10 19 Gauss) 
are predicted to exist Furthermore, it has been 

suggested that extremely strong magnetic fields (B ~ 
10 23 Gauss) are generated via the cosmological elec- 
troweak phase transition [7] and quark-hadron phase 
transition [Sj. Thus, QCD matter in strong magnetic 
fields with \eB\ ~ Aq CD is common in our universe. 

The analysis of chiral effective models has shown that 
under an external magnetic field, the chiral symmetry is 
always broken due to dimensional reduction. This phe¬ 
nomenon is called magnetic catalysis [mu and the chi¬ 
ral restoration temperature increases with the strength 
of the magnetic field. By contrast, the chiral restoration 
density decreases with a fixed low temperature fl6l - fl8} . 

A recent lattice simulation |jj| of QCD with an exter¬ 
nal magnetic field indicated that the critical restoration 
temperature decreases with a stronger magnetic field, 
which is called inverse magnetic catalysis or magnetic 
inhibition [joj - fjjl . The analysis presented in the present 
study does not provide an explanation of this effect, but 
instead it might suggest that the physics responsible for 
the magnetic inhibition at high temperature is beyond 
the physics captured by the Nambu-Jona-Lasinio (NJL) 
model. 

In this study, we analyze the NJL model, which de¬ 
scribes dynamical chiral symmetry breaking (DySB) [25l — 
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[27j |. in a strong magnetic field at finite temperature and 
finite density using the functional renormalization group 
(FRG) [28l - l33l | . The basic properties of this system such 
as the renormalization group (RG) flow equations and 
the fixed point structure have been investigated in many 
studies 0, Iddl - ifiOl ; see also the review paper [5lJ. We 
investigate the behavior of the RG flow of the four-Fermi 
coupling constant and we provide a detailed analysis of 
the phase diagram for the NJL model with an external 
magnetic field. 

The remainder of this paper is organized as follows. 
The formulation is given in the next section. Our results 
for structures of fixed points, the phase diagram, and the 
large -N non-leading effects are presented in Section Hill 
We summarize and discuss our results in Section HV1 


II. NJL MODEL AND ITS RG EQUATIONS 

In this section, we briefly introduce the RG flow equa¬ 
tion of the four-Fermi coupling constant and its struc¬ 
ture in an external magnetic field at finite temperature 
and density. The energy dispersion relation of a mass¬ 
less quark with electric charge q in an external magnetic 
field B = (0,0, B) is E 2 = (2n + l + s)\qB\ + pi, where 
s = ±1 is the Zeeman splitting due to the interaction be¬ 
tween the spin of a quark and B. The quantum number 
n = 0,1,2, • • • is the Landau level and in the particular 
case where n = 0 and s = — 1, this is called the Lowest 
Landau level (LLL). In the case of LLL, the dynamics of 
the quarks can be effectively described as a 1 + 1 dimen¬ 
sional system. 

We employ the following truncated effective action 
with Ul( 1) x Ur( 1) chiral symmetry in the Euclidean 
space, 


Ta = J dr J d 3 x \ip(0 + p, + qfi)ip 

Ga 


—+(#Y5V0 } 


(1) 


where the external vector potential A M is defined to give 
B = (0,0, B ) = rot A. The external field A M has several 
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representations due to the gauge degrees of freedom, e.g., 
the symmetric gauge: A M = (0,By/2,—Bx/2 1 0), or the 
Landau gauge: A M = (0, 0, Bx, 0). 

Theeffective action is governed by the Wetterich equa¬ 
tion [52L [Vi! . which in our case reads, 


d t T A [ip,ip] = -Tr 


d t RA 


-.(i.i) 


i?A 


( 2 ) 


where denotes the i-th (j-th) left (right)-hand side 
derivative of the effective action Ta with respect to ip (ip). 
The cut-off profile function Ra controls the shell momen¬ 
tum integration, thereby realizing the coarse-graining. 
We employ the 3d optimized cut-off function J5J, [55] , 


A-d 


Ra(p) = ii> - lj 0(1 - \p\/A). (3) 

The momentum mode integral at finite temperature in 
the external magnetic field takes the following form, 


d 4 p 


-A T 


oo 

E 

m =—oo 


\qB\ ^ f dp z 


27T 


1=0 


( 4 ) 


where m is the Matsubara mode number, and the factor 
2 on the left-hand side and cq = 2 — Si y o on the right-hand 
side are the spin-degeneracy factors. We also rewrite the 
Landau level as 2n + 1 s = 21 with l = 0,1, 2, • • •. 

Details of the derivation of the RG flow equations in 
this system were described previously [l8| . The RG equa¬ 
tions for the effective action m are reduced as follows: 


d t g = -2 g 
d t T = f, 
dtp = p, 
d t B = 2 B, 


g 2 (lN c J 0 (T,p,B)- Ji(T,p,B)) 


( 5 ) 

( 6 ) 

( 7 ) 

( 8 ) 


where g is the dimensionless rescaled four-Fermi coupling 
constant GaA 2 /27t 2 , T, p, and B are dimensionless exter¬ 
nal parameters, and dt denotes the derivative with re¬ 
spect to the dimensionless cut-off scale t = log(Ao/A). 
The threshold functions Jq and J\ are defined by 

I—M 

I Rl 2|?B| y- 

J 0 (T,p,B)='-^- J2 cnsJl^2l\qB\ 

1=0 

x {1 - n+ - n- - d t (n+ + n_)} , (9) 

I —WI 

I Rl 2|,B| /- 

J 1 (T,p,B)= 1 -^ J2 a^l-2l\qB\ 

x {iTTw{l- n *) + ojw{l- n -) 

- -——^d t n+ - —Edtn-l. (10) 
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FIG. 1: Dependence of the position of the UVFP g* on the ex¬ 
ternal magnetic field at finite temperature with the vanishing 
density in the large-IV leading approximation. 


The Gauss symbol |_£j denotes the greatest integer that 
is less than or equal to x. The Fermi-Dirac distribution 
functions n± are defined as 

7l± ~ e^( A± G + 1 _ e /3(!±M) + 1' ( ' 11 ^ 

The factor 4 and N c on the right-hand side of Eq. (0 de¬ 
note the number of degrees of freedom for the spinor and 
the color of the fermionic fields, respectively. The large- 
N leading approximation neglects the quantum correc¬ 
tions corresponding to the term J\ in Eq. ©. 


III. RESULTS 
A. Fixed Point Analysis 

In this subsection, we analyze the structure of the ul¬ 
traviolet fixed point (UVFP) g*, which satisfies fd(g*) = 0 
for the P function of g in Eq. The four-Fermi cou¬ 
pling constant g corresponds to the chiral susceptibility 
((ipip) 2 ). Therefore, the RG flow equation of g with the 
initial value go > g* diverges at a critical scale, where 
DySB with the second order phase transition turns on. 

In the case of strong B where the LLL approximation 
is effective with vanishing T and /i, the beta function of 
g is given by 


p g =-2g + N c g 2 \qB\, (12) 

which implies that the UVFP and the Gaussian fixed 
point (g* = 0) will become degenerate at large values of 
t. This property is due to the dimensional reduction [l8[ . 
In fact, the RG flow of g always diverges when go > 0. 
The thermal and density effects resolve the degeneracy 
of the fixed points. We investigate the dependence of 
the position of the UVPF on the temperature, chemical 
potential, and external magnetic field. The RG equa- 
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FIG. 2: Dependence of the position of the UVFP g* on the external magnetic field and finite density with fixed temperature 
(T/Ao = 0.0001) in the large -N leading approximation. The diagram on the right-hand side shows an enlargement of the small 
B region. 



FIG. 3: Dependence of the position of the UVFP g * of 
Eq. m on the external magnetic field with fixed chemical po¬ 
tential (p/Ao = 0.10) and several temperatures in the large -N 
leading approximation. 


T/Ao=0.0001 n/A«=0.1 



FIG. 4: Dependence of the position of the UVFP g* on various 
numbers for the color N c . The scale of the vertical axis is 
arbitrary. 


tion for the four-Fermi coupling constant in the large -N 
leading approximation is given by 

d t g = -2g + m c g 2 J 0 (f,jl,B), (13) 

where we set the following values: N c = 3 and q = (2/3)e. 
The UVFP g* of this equation is given by 


6J 0 (T,g,B)' 

The position of the UVFP with a finite temperature, fi¬ 
nite external magnetic field, and vanishing chemical po¬ 
tential is shown in Fig. [T] We can see that g* decreases 
monotonically as B increases for any temperature, and 
thus magnetic catalysis occurs. This result agrees with 
previous studies of chiral effective models (e.g., 0 )- 
The change in the position of the UVFP caused by 
finite density and a finite external magnetic field with 
a fixed low temperature (T/Aq = 0.0001) is shown in 


Fig. [2j For the smaller B region, the UVFP shown in the 
right-hand side panel of Fig. [2] oscillates. This behavior 
can be interpreted as the de Haas-van Alphen (dHvA) 
effect, which was observed in related studies, e.g., in color 
superconducting matter f5(|, holographic matter fl6j.[T7|. 
and in the NJL model with the mean field approxima¬ 
tion [57}. Thus, the oscillatory behavior is derived from 
the processes at each Landau level 2\qB\n that cross the 
Fermi surface /i 2 . The end point of the oscillation corre¬ 
sponds to the case where the first Landau level (n = 1) 
just overlaps with the Fermi surface. Therefore, only 
the LLL can make a contribution beyond the oscillatory 
region. For example, in case where p./Ao = 0.2 the oscil¬ 
lation ceases at eB/ A§ ~ 0.03. These values satisfy the 
relationship 2\qB\ = 2||eB| ~ /i 2 . The dHvA effect dis¬ 
appears at higher temperature (see Fig. The Landau 
level and the Fermi surface do not depend on the num¬ 
bers of the flavor and color, and thus this phenomenon 
is not affected by them, as shown in Fig. 0] 



































































































































4 


H /4=0.25,T/4=0.06 n /4=0.25,T/4=0.06 




FIG. 5: Comparison of the position of the UVFP g* in the large-IV leading case and non-leading case with N c = 1 (left) and 
A c = 2 (right). 


large-N leading 



|i/Ao 


large-N non-leading 



jn/Ao 


FIG. 6: Chiral phase diagram of the fi-T plane. The figures on the left and right are the large-IV leading case and non-leading 
case, respectively. 


After the oscillation ceases, the UVFP g* increases 
with B , which means that inverse magnetic catalysis oc¬ 
curs. The inverse magnetic catalysis at low temperature 
is reproduced qualitatively by analyzing the fixed point 
structure jl6l - [l8| . This phenomenon cannot be observed 
in a lattice simulation at present because it is affected by 
the so-called sign problem. 

Moreover, we note that the UVFP decreases for larger 
values of A. It appears that the effect of finite density de¬ 
couples with a larger external magnetic field because the 
magnetic field at the peak is proportional to the chemi¬ 
cal potential. At present, there is no clear explanation of 
this behavior. 

We investigate the impact of the large-A non-leading 
corrections using Eq. ©, where the term J\ is included. 
A comparison of the position of the UVFP g* in the large- 
A leading case and the non-leading case with A c = 1 
and A c = 2 is shown in Fig. [5] The large-A non-leading 
corrections make the UVFP g* larger in the whole region 
of B. However, the qualitative magnetic behavior of the 


UVFP does not change greatly after the inclusion of the 
large-A non-leading effects. 


B. Phase diagram 

The phase diagrams for the fi — T plane with fixed ex¬ 
ternal magnetic field are shown in Fig. [6] We evaluate the 
RG equation of the inverse four-Fermi coupling constant 
in order to obtain the four-Fermi coupling constant at 
the infrared limit [5(j. Note that the phase boundaries in 
these phase diagrams do not necessarily indicate the true 
chiral phase transition point. Instead, these boundaries 
represent the instability of the chiral symmetric vacuum, 
which is the phase transition point only for the second 
order transition. For the first order phase transition, the 
true transition point moves to the symmetric side. We 
set the dimensionless rescaled inverse four-Fermi coupling 
constant to 1/g = 0.182 at the initial scale t = 0, for 
which the RG flow at vanishing temperature and density 
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shifts to the broken phase at the infrared scale. 

First, we discuss the large -N leading case. In the high 
temperature region, the phase boundaries move to the 
higher temperature side as B increases, i.e., magnetic 
catalysis occurs. By contrast, inverse magnetic catal¬ 
ysis occurs in the low temperature region. The phase 
boundaries then move toward the lower density side as 
B increases. Next, for larger values of B 1 they turn back 
toward the higher density side. This behavior is also ob¬ 
served in the mean-field approximation [57| . It appears 
that this behavior is due to competing effects between 
magnetic catalysis and inverse magnetic catalysis. 

Next, we discuss the large-iV non-leading effects. At 
low temperature, the phase boundary turns toward the 
lower density side, even with the vanishing B. This oc¬ 
curs due to the singularity at the Fermi surface (see ]50}). 
At vanishing temperature and finite density, the non¬ 
leading correction becomes singular at /j, = A, so the 
non-leading effect becomes larger than the leading one in 
the low temperature region. As shown by Eq. (j5j , the 
non-leading term has a negative sign in the RG equation 
of the four-Fermi coupling constant, i.e., the non-leading 
effect makes the phase more symmetric. At finite eB , we 
still see that the phase boundaries move back and forth 
as the external magnetic field increases in the low tem¬ 
perature and high density region. Thus, the non-leading 
effects do not change the behavior attributable to the 
external magnetic field. 


IV. SUMMARY AND DISCUSSION 

In this study, we investigated the dependence of the 
UVFP on the four-Fermi coupling constant in the NJL 
model at finite temperature and density under an ex¬ 
ternal magnetic field by using the FRG. The UVFP de¬ 
creases monotonically as the magnetic field increases at 
finite temperature and vanishing chemical potential, and 
thus magnetic catalysis occurs. At finite chemical po¬ 
tential and a fixed low temperature, the UVFP oscil¬ 
lates depending on the external magnetic field due to the 
dHvA effect. This effect vanishes at higher temperatures. 
The UVFP increase as the magnetic field increases af¬ 
ter the oscillatory region, which means that inverse mag¬ 
netic catalysis occurs. However, an even larger external 
magnetic field changes the inverse magnetic catalysis into 
magnetic catalysis. 

We also investigated the chiral phase diagram. Mag¬ 
netic catalysis is observed at high temperature and low 
density. However, at low temperature and high density, 
inverse magnetic catalysis occurs with a large external 
magnetic field. For a larger external magnetic field, the 
phase boundaries move back to the large density side, and 
thus magnetic catalysis occurs. The large-TV non-leading 
effects do not change the qualitative behavior of our sys¬ 
tem dramatically, although the phase boundary moves 
toward the symmetric side at low temperature and high 
density in eB = 0 due to the singularity at the Fermi 


surface (Hoj . 

In order to investigate the cut-off scheme dependence, 
we analyzed the system with the Id optimized cut-off 
function {58} , 


Ra(Pz) = Wzlz 



l)fl(l-|p*|/A). 


(15) 


The threshold functions are shown in Appendix [A] We 
found that the qualitative behaviors do not change, such 
as the dependence of the UVFP on thermal effects and 
the shapes of the phase boundaries. Clearly, the values of 
the critical temperature, density, and external magnetic 
field change because the NJL model is itself an unrenor- 
malizable theory. It appears that the behaviors of the 
UVFP and the phase diagram determined in this study 
are quite stable relative to the cut-off profile. 


We also comment on the momentum-dependent cou¬ 
pling constant, i.e., the non-local vertex G(p). The non¬ 
local vertex is partly included through the large-TV non¬ 
leading diagrams, and thus they can be considered to 
represent some of the fluctuations in mesons (18]. The 
large -N leading effect is much larger than the non-leading 
effects, so our results only change slightly. After includ¬ 
ing the momentum-dependence of the four-Fermi cou¬ 
pling constant, the phase boundary is expected to move 
toward lower temperature and density because the chiral 
symmetry tends to be restored by mesonic fluctuations. 


We hope that these analyses motivate more elab¬ 
orate studies in the future using the re-bosonization 
method [59l - (63[ or the weak solution method 64}. In¬ 
deed, more precise analyses should be performed. 
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Appendix A: Threshold function with the Id 
optimized cut-off function 


We give the threshold functions using the Id optimized 
cut-off function (fT5l) in the RG equation of the four-Fermi 
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coupling constant. We have 

WB\ 


I 0 (TJt,B) = 


-1 


C(3 

V2\eB\ 3 / 2 


~J 2 ai 

1=0 

h(T,fi,B) = 


(n + +n_) <9 t (n + + n_) 


m 


E 

;=o 


ai 


Wl? \2 


- n- - 


Ul (^) 2 V 2 

5 t ?z + <9 t n_ 


o ~ n + 


(Al) 


(A2) 


where (^(s,x) is the Hurwitz zeta function, n± are given 
in Eq. (flTT) . and = e)±/2 with ef = l + 2l\qB\. The RG 
equation of the four-Fermi coupling constant is obtained 
by replacing Jq and J\ with /q and I±, respectively. 
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